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ABSTRACT 

It  Is  observed  that  normally,  during  the  polar  night,  24- 
hour  periods  with  a  clear  sky  and  light  to  moderate  winds  have 
their  lowest  temperatures  in  the  middle  of  the  period,  and  the 
highest  temperatures  occur  at  the  beginning  and  end.   For 
cloudy  and  windy  periods,  the  shape  of  the  temperature  curve 
is  reversed,  with  the  maximum  temperature  occurring  near  the 
middle  of  the  period.   Despite  many  efforts  by  various  invest- 
igators no  physical  explanation  could  be  found  for  this  vari- 
ation. 

Hisdal  discovered  that  for  clear  skies  and  low  wind  speed, 
the  same  shape  of  temperature  variation  curve  occurs,  whether 
the  starting  time  is  0000  hr,  1200  hr,  or  1400  hr.   Thus,  the 
abnormal  behavior  of  the  temperature  variation  appears  to  be 
influenced  primarily  by  the  statistical  bias  involved  in  the 
selection  criteria. 

This  investigation,  using  six  years  of  surface  data,  di- 
vides the  diurnal  temperature  variation  for  various  special 
types  of  days  into  two  parts,  the  meteorological  or  diurnal 
variation  and  the  superimposed  statistical  bias.   Finally,  the 
meteorological  diurnal  temperature  variation,  based  upon  all 
observations  during  the  six  year  period,  is  determined. 

The  writer  wishes  to  express  his  appreciation  for  the 
advice  and  guidance  given  him  by  Professor  W.  van  der  Bijl  of 
the  U.  S.  Naval  Postgraduate  School  in  this  investigation. 
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1 .   Introduction 

In  his  meteorological  investigations  with  the  Norwegian- 
British  -  Swedish  Antarctic  expedition  of  1949  -  1952,  V. 
Hisdal  [l  j  discusses  a  much  pondered  peculiarity  of  the  diur- 
nal temperature  variation  during  the  polar  night.   For  the 
various  polar  stations  investigated^  days  with  clear  skies 
and/or  light  to  moderate  winds  had  their  minimum  temperatures 
near  1200  GCT,  whereas  for  the  cases  of  overcast  skies  and/or 
windy  days  the  opposite  type  of  diurnal  variation  occurred; 
that  is9  a  maximum  temperature  near  1200  GOT, 

Concerning  this  strange  temperature  behavior^  Simpson 

[2,  p. 72]  wrote  in  1919s 

I  have  spent  more  time  over  this  problem  than  over 
any  other  discussed  in  this  work  ,,,»«>  the  investi- 
gation has  led  only  to  negative  results , 

In  1951  Court  [3]  could  still  not  explain  this  tempera- 
ture behavior.   He  wrote % 

Why  should  the  maximum  temperature  on  clear  days 
during  the  antarctic  night  occur*,  almost  without 
exception9  in  the  hour  after  midnight?  Simpson 
pondered  this  problem  longer  than  any  other  of 
the  many  he  investigated  in  compiling  the  meteor- 
ological reports  of  the  last  Scott  expedition;, 
without  developing  a  satisfactory  hypothesise 

Rouch  C4.J  also  found  a  dally  temperature  variation  during 

the  dark  seasony  based  on  an  analysis  of  the  data  both  from 

the  Charcot  expeditions  and  from  the  Bird  expeditions.   He 

wrote  that  the  average  daily  variation 

is  a  simple  curve  with  a  maximum  around  1  or  2 
A.M.  and  minimum  after  noon  eeee  I  cannot  give 

a  satisfactory  explanation. 


Hisdal's  explanation  for  this  phenomenon  was  simply  that 
a  physical  explanation  was  not  necessary  at  all.   It  could  be 
explained  in  terms  of  an  "automatic"  daily  variation  which  re- 
sults from  the  nature  of  the  type  of  day  selected.  For  ex- 
ample, Hisdal  [5 9   p. 104)  has  shown  that  if  a  type  of  day  is 
characterized  by  the  presence  of  a  temperature  inversion,  it 
will  usually  have  a  lower  average  temperature  than  the  non- 
inversion  days.   This  means  that  if  the  average  diurnal  tem- 
perature curve  is  computed  using  only  the  inversion-type  days 
(having  lower  temperatures) ,  the  non-inversion  days  (having 
higher  temperatures)  will  automatically  tend  to  give  a  contri- 
bution near  both  edges  of  the  24-hour  time  span.  The  result 
of  this  automatic  daily  variation,  shown  schematically  in  fig- 
ure 1 ,  is  that  it  superimposes  a  statistical  bias  upon  any 
given  diurnal  variation  that  might  be  present.  Moreover,  the 
magnitude  of  the  bias  could  be  large  enough  to  mask  the  diurnal 
variation  completely. 

An  interesting  fact  is  that  Hisdal's  explanation  for  the 
peculiar  behavior  of  the  diurnal  temperature  variation  is  ap- 
parently an  independent  discovery  of  the  theory  that  Bartels 
C6H  described  in  his  Gottingen  dissertation  of  1922. 

Hisdal  demonstrated  this  statistical  bias  explicitly  by 
considering  two  diurnal  temperature  variation  curves  for 
polar  stations 9    based  on  data  for  dark  season  days  with  wind 
speeds  less  than  15  miles  per  hour  and  cloudiness  less  than 
five- tenths.   The  first  curve  begins  at  0000  GOT  and  has  a 
minimum  temperature  at  1200  GOT.   The  second  begins  at 
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Figure  1 
Schematic  diagram  showing  the  statistical  bias  effect 


1200  GOT  and  has  a  minimum  at  approximately  0000  GOT.   If 
there  were  a  uniform  diurnal  temperature  variation,  the  min- 
imum temperatures  should  always  occur  at  the  same  local  time. 
Moreover,  varying  the  initial  time  of  the  graph  should  not 
change  the  time  of  occurrence  of  the  minimum.  A  logical  ex- 
planation of  the  difference  in  these  diurnal  temperature  vari- 
ations is  that  the  minima  were  primarily  due  to  the  statistical 
bias>  not  by  diurnally  different  physical  processes. 

W.  van  der  Bijl  [7,  p.  149]  also  supports  the  theory  of  a 
"statistical  bias"  in  his  investigation  of  the  diurnal  pressure 
variation  in  northwest  Europe.  His  study  discloses  a  similar 
automatic  diurnal  variation  due  to  data-selection  bias.   The 
basis  of  selection  was  that  each  day  had  to  be  in  the  month  of 
December  and  have  at  least  5.7  hours  of  insolation.  These  cri- 
teria resulted  in  a  diurnal  trend  with  low  pressure  at  the 
beginning  and  end  of  the  period,  with  high  pressure  in  the 
middle. 

Thus,  it  has  been  shown  that  the  existence  of  a  "statis- 
tical bias"  or  "automatic  diurnal  variation"  gives  a  reason- 
able explanation  of  the  behavior  of  the  temperature  trend  at 
polar  stations  during  the  dark  season.  However,  Hlsdal  in 
his  investigation  did  not  deduce  the  mean  meteorological  diur- 
nal variation  of  temperature,  even  though  a  knowledge  of  the 
effect  of  the  statistical  bias  made  this  deduction  possible. 

The  purpose  of  this  study  is  to  compute  the  mean  meteor- 
ological or  diurnal  temperature  variation  during  the  polar 
night.  As  a  secondary  purpose,  the  temperature  variations 


induced  by  the  statistical  case-selection  procedure  will  also 
be  computed. 


2.   Data  and  Processing 

The  data  used  in  this  investigation  were  obtained  from 
the  National  Weather  Records  Center,  Asheville,  North  Carolina. 
The  meteorological  information  consisted  of  punched  IBM  cards 
of  surface  data  at  three-hourly  intervals  from  McMurdo  Sound, 
Antarctica  (78S,  166E),  for  the  years  1956  through  1961.  Fig- 
ure 2  shows  the  geographical  location  of  the  area  concerned. 

The  duration  of  the  polar  night  at  McMurdo  sound  was  de- 
termined by  the  use  of  the  solar  altitude  curves  prepared  by 
LeVasteau  [8j  •   These  indicated  that  the  duration  of  the  dark 
season  was  from  20  March  through  20  September.   However,  this 
period  does  not  eliminate  twilight  effects  and  incoming  solar 
radiation  due  to  unusual  atmospheric  refraction  conditions. 
Therefore,  the  shorter  period,  15  May  through  15  August,  was 
used  in  order  to  insure  complete  darkness. 

In  processing  the  punched  data  cards ,  the  IBM  82  card 
sorter  was  used  to  select  observations  for  the  15  May  through 
15  August  period.   Then  the  card  data  was  transferred  to  a 
magnetic  tape  by  use  of  the  IBM  1401  printer.   The  magnetic 
tape  was  subsequently  used  in  conjunction  with  a  program  for 
the  Control  Data  Corporation  1604  digital  computer,  in  order 
to  determine; 

(1)  the  mean  diurnal  variations  of  temper- 
ature, cloudiness,  atmospheric  pressure, 
and  wind  direction  for  various  years; 

(2)  the  same  as  (1),  except  for  various 
starting  times  of  the  period; 


(3)  the  same  as  (2),  except  for  periods 
other  than  24  hours; 

(4)  the  total  six-yearly  mean  variations 
as  stated  in  (1),  (2),  and  (3)  above; 

(5)  and  the  Lamont  correction  for  the  above 
mentioned  curves. 


Figure  2 
Geographical  location  of  McMurdo  Sound,  Antarctica 
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3.   The  Lamont  correction 

The  Lamont  correction  was  mentioned  in  section  2  as  one 
of  the  computer  subroutines.   The  purpose  of  this  correction 
is  to  make  the  starting  value  of  a  continuous  parameter  and 
the  final  value  of  the  parameter  equal  to  each  other.   This 
technique  simplifies  the  study  of  any  cyclical  variations  and 
is  especially  valuable  in  determining  a  diurnal  variation. 

A  more  complete  discussion  of  the  Lamont  correction  is 
presented  by  Conrad  and  Pollak  [9]    .   For  this  study  the  cor- 
rected temperature  was  computed,  using  equation  (2) 1 

TiL  =  3;  +  Ti       (2) 

31  =  '(t°-Il)     (3) 
n 

where   T;L  =  corrected  temperature 

T;  =  sequence  of  three-hourly  temperature  observations 

g.  =  Lamont  correction. 
A  typical  application  of  the  Lamont  correction  follows. 
For  this  purpose,  use  is  made  of  equations  (2)  and  (3),  and 
the  hypothetical  temperature  observations  given  in  table  1. 
The  computed  values  of  T-iL    are  summarized  in  the  lower  half  of 
table  1 • 


given 

To    - 

10 

T4  = 

2 

n     = 

4 

thus 

3;  = 

(H^1    =("2)i> 

and 

t;l  = 

(-2);  +  t,  . 
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Table  1 

Some 
the 

hypothetical  temperature  observations (Tj  )  and 
resulting  Lamont-corrected  temperatures (Til  ) 

Hour 

of  observation  (GOT)    00 

06 

12 

18   24 

i 

0 

1 

2 

3    4 

*i 

2 

3 

6 

9   10 

Si 

0 

-2 

-4 

-6   -8 

TiL 

2 

1 

2 

3    2 

It  must  be  noted  that  the  temperatures  as  corrected  by 
the  Lamont  procedure  are  quite  different  than  the  original, 
uncorrected  temperatures.   However,  since  the  present  investi- 
gation deals  only  with  the  meteorological  temperature  cycle 
(i.e.  the  T;j_  values),  actual  values  of  temperature  are  not  im- 
portant. 

If  the  actual  temperature  values  were  desired,  it  would 
be  necessary  to  maintain  a  table  of  corrections,  as  well  as 
the  corrected  values. 


4,   Theoretical  discussion  of  the  method  used  to 
separate  the  meteorological  diurnal  variation 
and  the  statistical  variation 

As  previously  explained,  the  apparent  diurnal  temperature 
variation  observed  during  the  polar  night  for  special  types  of 
days  may  be  considered  to  be  the  sum  of  the  meteorological  var- 
iation and  an  effect  resulting  from  statistical  bias.   More- 
over, the  major  purpose  of  this  study  is  to  determine  an  esti- 
mate of  these  two  components. 

In  order  to  illustrate  how  this  dichotomy  may  be  performed, 
it  will  be  helpful  to  consider  first  how  a  particular  observed 
daily  temperature  variation  is  "built  up"  by  combining  a 
meteorological  influence  with  a  statistical  one.  After  the 
total  diurnal  variation  has  been  constructed,  it  will  once  again 
be  separated  into  its  original  components.   This  example  will 
serve  to  illustrate  how  the  actual  data  for  this  study  may  be 
treated  to  derive  the  meteorological  and  statistical  variations. 

First,  consider  that  the  set  of  T-,|_  values  listed  in 
table  1  are  representative  of  the  meteorological  variation. 
These  temperatures  will  remain  fixed  for  the  particular  hours 
for  which  they  were  reported  (e.g.,  the  meteorological  influ- 
ence for  0600  GOT  will  always  be  1F,  regardless  of  the  starting 
time  of  the  period) . 

Next,  consider  a  purely  hypothetical  statistical  vari- 
ation as  listed  in  table  2.  As  previously  pointed  out  by 
Hisdal,  the  statistical  bias  is  not  fixed  at  a  particular 
value  for  a  particular  observation  time,  but  depends  upon  the 
starting  time  of  the  period  selected.   To  illustrate  this, 
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assume  first  that  the  starting  time  of  the  24-hour  period  in 
table  2  is  0000  GOT.   Then  the  statistical  bias  for  0600  GOT 
would  be  1P,  since  it  would  be  the  second  observation  reported. 
On  the  other  hand,  if  0600  GOT  were  the  starting  hour  of  the 
24-hour  period,  the  0600  GOT  bias  would  be  3F.   Thus,  the  sta- 
tistical influence  may  be  thought  of  as  a  wave  which  is  free 
to  move  through  the  diurnal  or  meteorological  wave  with  con- 
tributions dependent  upon  the  starting  time  of  the  24-hour 
period. 

Table  2 

Hypothetical  statistical  bias 
for  each  6-hourly  observation 
time  in  a  24-hour  period 

Observation  time  i   Temperature (F) 
1st  3 

2nd  1 

3rd  0 

4th  1 

5th  3 

Now,  using  the  meteorological  influence  from  table  1  and 
the  statistical  influence  from  table  2,  the  "total"  tempera- 
tures will  be  computed  for  various  starting  times,  as  shown 
in  table  3. 

The  values  of  total  temperature  for  the  various  24-hour 
periods  are  summarized  in  table  4.   The  smaller  numbers  adja- 
cent to  the  total  temperature  are  the  meteorological  and 
statistical  contributions  to  the  total,  as  follows. 
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Table  3 

Total  temperature  variation  resulting  from  the  statistical 
and  meteorological  variations  from  tables  1  and  2 
respectively j  for  various  starting  times 

■—— ——————      J     ■  —— — — — — — MJ— — — 

Starting  time:  0000  GOT 

Hour:  00   06    12  18   00 

Statistical:  3    10  13 

Meteorological:  2    12  3    2 

Total:  5    2    2  4    5 

Starting  time:  0600  GOT 
Hour:  06   12   18   00   06 

Statistical:  3  10  13 
Meteorological:  12  3  2  1 
Total:  4    3    3    3    4 

Starting  time:  1200  GOT 

Hour:             12   18   00  06   12 

Statistical:       3    10  13 

Meteorological:    2    3    2  12 

Total:            5    4    2  2    5 

Starting  time:  1800  GOT 

Hour:  18   00   06   12   18 

Statistical:  3    10    13 

Meteorological:  3    2    12    3 

Total:  6    3    13    6 
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g 

tin;  where  t  =  total 

s  =  statistical  influence 
m  =  meteorological  influence. 
The  total  temperature  variations  in  table  4  are  typical 
of  those  obtained  from  the  CDC  1604  computer  program.   They 
are  plotted  in  figure  3  in  order  to  show  the  effect  the  stat- 
istical influence  has  in  shifting  the  time  of  the  daily  min- 
imum for  various  starting  times  of  the  hypothetical  24-hour 
period. 

Now,  the  next  step  in  this  discussion  is  to  perform  a 
dichotomy  upon  the  total  temperature  values  found  in  table  4 
in  order  to  determine  the  original  values  of  meteorological 
and  statistical  variation  contributions.  A  more  general  proof 
of  the  separation  method  is  given  in  section  5.   Specifically, 
the  method  to  be  used  is  the  same  as  that  utilized  by  van  der 
Bijl  in  his  investigation  of  diurnal  pressure  variations.   It 
consists  of  averaging  the  values  of  table  4  in  either  of  two 
ways:  by  columns,  to  determine  the  meteorological  influence, 
or  diagonally,  to  determine  the  statistical  influence.   Once 
either  of  these  averages  has  been  computed,  the  other  may  be 
determined  by  subtraction  from  the  total  values.   These  methods 
are  shown  below.   The  notation  used  is  as  follows:  tz  ^ 

t  =  temperature  (P) 

0  =  m;  refers  to  meteorological  influence 
$     -     s;  refers  to  statistical  influence 
hh  =  hour  of  observation  (GOT). 
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Table  4 

Hypothetical  total  temperatures 
for  various  24  hour  periods 


Hour  (GOT)   00   06   12   18   24/00  06   12   18 
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Figure  3:  Hypothetical  total  temperature 
variation  curves  for  various  starting 
times  of  the  period,  based  upon  the 
values  found  in  table  4. 
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Averaging  by  columns,  and  subtracting  the  constant  value 
of  the  statistical  bias,  results  in  the  meteorological  varia- 
tion. 

+  ( L±l±2z±£\-(  l±±±-Q-±±)  -ML 

tmfi0={~      3         j    {         4       j  -    -$    ^ 

=  Z 

+         -    /  H  4+Z+l\_  (  I  +  3  +  \  +  0\.  ±-  S 

=  1 

/  z*  3  +5  +  3)     (Q  +  I  +  3  +  t).    13-5 
.  /4+3  44tQ-//f  °*  /  *3)  =  IZ-i 

=  3 

f  /^3tZ^3)      /34    |^0+|\.I3,^ 

Averaging  diagonally  (assuming  the  24-hour  period  begins 
at  0000  GOT),  and  subtracting  the  constant  value  of  the  mete- 
orological influence,  results  in  the  statistical  influence. 

tSj00   =  (     4     J      \  +  J        T      ? 

=  3 


tM.-  (&^±ii;-(t^Hi^;-f- 


-? 
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=  0 

=  1 

Vf  =  (  — 4 y   I 3 J    t   ? 

=  3 

In  this  example  we  have  recovered  the  values  of  the  mete- 
orological and  statistical  influences,  assumed  originally  in 
constructing  the  values  of  total  temperature. 

One  difficulty  in  applying  these  methods  to  actual  data 
lies  in  the  fact  that  the  constant  terms  which  are  subtracted 
from  the  averages  of  total  temperature  are  not  immediately 
apparent.   This  difficulty,  however,  is  of  no  great  concern, 
when  the  temperature  trend  and  not  the  actual  temperature  per 
se  is  under  investigation. 

For  the  sake  of  completeness,  the  next  section  will 
illustrate  how  the  use  of  these  constants  may  be  eliminated. 
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5.    Elimination  of  the  constant  term,  when  computing 
the  statistical  bias 

The  purpose  of  this  section  is  to  illustrate  how  the  con- 
stant term  may  be  eliminated,  when  computing  the  statistical 
bias.   The  only  known  temperature  values  are  those  of  total 
temperature,  found  in  table  4,   The  hourly  meteorological  con- 
tributions are  unknown. 

The  notation  used  is  as  specified  in  section  4,  except 
for  the  following: 

Cm  =  a  constant  which  is  dependent  upon  the 

meteorological  influence, 
Cs  =  a  constant  which  is  dependent  upon  the 
statistical  bias. 
The  column  averages,  computed  in  section  4,  become: 
t    =  '3  -  r 

t^   -  q  - r 


Next,  if  the  averaging  is  done  diagonally,  using  OOOO 
GrCT  as  the  starting  time,  the  following  temperatures  result: 

%      =  ^  -  C 

w  s^ofe  -   4    L.  m  5  m     ^ —  j 
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Making  the  appropriate  substitution  for  Cm  in  the  last 
set  of  equations_, 

-  4"*  i\^  4  f  -v  v, — — — ^    *-j 

"5,00  -   ^     J6       it — ~ 

-'•  5tSj0o-tSj(>i>-%%-tSjti   s    7    ,  (4) 

Similarily,  for  the  starting  times  06,  12  and  18  h,  we 
obtain,  respectively: 

tf,«>+ is/**  -is,ix'3tSJS=   1  .  (7) 

Equations  (4)  through  (7)  comprise  a  set  of  four  linear, 
homogeneous  equations  in  four  unknowns,  which  can  be  readily 
solved  using  the  usual  mathematical  techniques.   Such  a  solu- 
tion produces  the  following: 
V  =  3 

H*   =  ° 

"^  5,2-4   =   *S,00  =   3  . 
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Thus,  it  is  seen  that  if  specific  values  of  statistical 
bias  or  meteorological  effect  are  desired,  a  set  of  homogeneous 
equations  must  be  solved.   However,  as  was  stated  in  section 
4,  this  present  study  does  not  require  specific  temperature 
values,  but  rather  their  relative  values.  Therefore,  the 
constants  0S  and  C^  will  not  be  removed  from  the  computed 
diurnal  averages. 
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6.   General  discussion  of  the  actual  results 

In  general,  the  results  of  this  study  can  be  separated 
into  four  parts.   The  first  is  the  set  of  24-hour  temperature 
variation  curves,  figures  4  through  25.   These  curves  were 
drawn  in  order  to  determine  to  what  degree  statistical  bias 
is  apparent,  when  dealing  with  actual  data.  Also,  in  these 
curves  the  criteria  used  in  selection  of  the  type  of  day  to 
be  investigated  were  varied.   The  extreme  cases  considered 
were  (1)  wind  speed  greater  than  15  kt>  cloudiness  greater 
than  five-tenths;  and  (2)  wind  speed  calm,  cloudiness  equal  to 
zero. 

The  second  set  of  curves,  figures  26  through  30  are  the 
48-hour  variations  for  various  pairs  of  days  selected.   The 
selection  criteria  are  comparable  to  those  specified  for  the 
24-hour  variation  curves.   The  major  purpose  of  the  48-hour 
curves  is  to  show  that  the  statistical  bias  effect  is  not  res- 
tricted to  a  24-hour  period. 

The  third  set,  figures  31  through  34,,  comprise  the  18-hour 
and  the  72-hour  variations.   These  curves  are  also  included  in 
the  study  to  show  how  periods  of  various  length  effect  the 
statistical  bias. 

The  fourth  set  of  curves,  figures  35  through  42,  show  the 
relationship  between  temperature  variation  and  variations  of 
atmospheric  pressure,  cloudiness,  and  wind  direction.   The  pur- 
pose of  this  set  is  to  aid  in  giving  some  explanation  for  the 
peculiarities  found  in  the  daily  temperature  trend. 

In  all  of  these  curves  the  abscissa  is  the  hour  of  obser- 
vation in  Greenwich  Civil  Time,  which  is  the  same  time  scale 
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that  Hisdal  used  in  his  investigation.   It  is  convenient  to 
use  Greenwich  Civil  Time  rather  than  local  time,  since  the 
reporting  times  of  surface  observations  are  in  GOT. 

If  it  is  desired  to  obtain  the  diurnal  variation,  based 
on  local  time,  add  11  hours  to  the  observation  times  in  Green- 
wich Civil  Time.   However,  it  must  be  remembered  that  local 
time  loses  its  significance  during  the  polar  winter  season, 
since  the  noon  hour  has  no  more  sunshine  than  the  "witching 
hour" • 
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7.    The  diurnal  temperature  variation  for  special 
s e lections  of  days 


Table  5  summarizes  the  comparison  between  the  expected 
curvatures  of  the  diurnal  variations  for  selected  types  of 
day  and  various  starting  times  of  the  period  considered.  The 
expected  curvature  is  based  on  the  discussion  in  section  1. 
There  it  was  pointed  out  that  days  with  high  wind  speed 
and/or  cloudy  skies  will  have  above-normal  average  values  of 
temperature,  resulting  in  a  convex  diurnal  temperature  curve. 
Days  with  low  wind  speed  and/or  clear  skies  have  below=normal 
mean  temperatures,  resulting  in  a  concave  diurnal  temperature 
curve.   In  table  5  the  shape  of  the  daily  meteorological 
course  of  temperature  agrees  favorably  with  the  expected  shape. 
However,  as  the  criteria  are  made  more  restrictive,  the  agree- 
ment  becomes  worse ,  due  to  the  increasing  contributions 
arising  from  random  effects  as  the  sample  size  is  reduced. 
(The  effects  of  randomness  are  considered  in  section  11.) 

Table  6  contains  a  summary  of  the  expected  shape  and 
actual  shape  of  the  statistical  bias  curves  for  the  various 
types  of  day  selected.   The  table  also  contains  information 
concerning  the  time  of  maxima  or  minima  in  the  mean  meteoro° 
logical  curves  and  the  range  of  temperature  fluctuations 
involved.   It  is  seen  that  the  actual  shape  of  the  statistical 
bias  curves  agrees  with  the  expected  shape  quite  well.   How~ 
ever,  as  in  table  5S   making  the  selection  criteria  more 
restrictive  produces  more  pronounced  random  errors s    due  to 
the  reduction  of  the  sample  size.   For  instances  figure  21 
is  based  upon  only  seven  days.   The  result  is  that  no 
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Table  5 

Summarized  results  of  the  total  temperature 
variation  curves  for  various  starting;  times 
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organized  statistical  influence  Is  apparent  in  this  case.   In 
contrast,  figure  7  is  based  on  174  days  and  the  statistical 
bias  has  a  definitely  concave-shaped  curve. 

One  peculiarity  that  occurs  in  the  set  of  temperature 
variation  curves  is  the  2100  GOT  maximums  appearing  in  the 
statistical  bias  curves.   This  maximum,  which  would  occur  at 
0800  local  time,  gives  an  unnatural  appearance  to  the  bias 
curves.   Therefore ,  considerable  time  was  spent  in  determining 
the  cause  of  the  peculiarity.   One  resulting  conclusion  is 
that  the  2100  GOT  maximum  is  associated  with  the  criterion  of 
calm  winds,  since  it  appears  in  all  curves  in  which  calm  winds 
are  included,  except  those  of  figure  21 .   Here  the  maximum 
occurs  at  1800  GOT.   However,  this  latter  figure  is  based  on 
only  seven  days  and  one  should  be  careful  not  to  attribute 
undue  significance  to  such  a  small  sample. 

An  attempt  was  made  to  exclude  the  calm  days  in  computing 
the  mean  diurnal  temperature  variation  in  order  to  determine 
whether  the  2100  GOT  maximum  would  disappear  from  the  statisti- 
cal bias.   The  results  are  contained  in  figures  5  and  23 $    in 
which  the  maximum  does  not  appear.   However,  the  criterion  for 
figure  23  (wind  speed  between  two  and  15  kt)  was  too  restrictive 
and  only  four  days  could  be  used. 

Further  investigation  revealed  how  the  calm-day  observa~ 
tions  affected  the  statistical  bias  curves.   This  matter  is 
explained  in  more  detail  in  section  11.   It  will  suffice  the 
present  discussion  to  say  that  the  maximum  in  the  statistical 
bias  curve  resulted  from  calculating  the  bias  through  the  use 


26 


of  a  total  temperature  curve  which  had  an  above-normal  random 
contribution  at  2100  GOT. 

It  is  interesting  to  note  that  the  1 8-hour  statistical 
bias  curve,  figure  32,  does  not  contain  a  maximum  comparable 
the  the  2100  GOT  maximum.   Therefore,  it  is  doubtful  that  such 
a  maximum  would  occur  with  an  independent  selection  of  data. 
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Figure  4:  Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speed  ^15  let;  cloudiness  ^5/10) 
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Figure  5:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  (Met.)  and 
the  statistical  bias  variations. 

(wind  speed  ^15  let;  cloudiness 2*5/10) 
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Figure  6:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speed  £15  kt;  cloudiness-^  5/10) 
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Figure  7:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  (Met.)  and 
the  statistical  bias  variations. 

(wind  speed ^15  kt;  cloudiness  £  5/10) 
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Figure  8:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speeds  8  kt;  cloudiness^  5/10) 
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Figure  9'      Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  the 
statistical  bias  variations. 

(wind  speeds  8  let;  cloudiness£5/lO) 
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Figure  10:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speed  ^5  let;  cloudiness^5/l0) 
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Figure  11:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  the 
statistical  bias  variations. 

(wind  speed £5  ktj  cloudiness^5/lO) 
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Figure  12:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period. 
(wind  speeds  4  kt;  cloudiness^  5/10) 
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Figure  13:  Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  statis- 
tical variations. 

(wind  speed  ^4  let;  cloudiness^5/l0) 
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Figure  14:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speeds  3  let;  cloudiness^ 5/10) 
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Figure  15:  Comparison  of  the  total  temperature 
variation  with  the  meteorological  and 
statistical  bias  variations. 

(wind  speed  &J>   let;  cloudiness  <{■  5/10) 


-27.0 


09 

HOUR 


12    15 

OF  OBSERVATION 
(GOT) 


39 


Figure  16:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speed -i  2  kt;  cloudiness^5/l0) 
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Figure  17:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  the 
statistical  bias  variations. 

(wind  speeds  2  kt;  cloud  iness-^5/10) 
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Figure  18:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speed  =  0  kt;  cloudiness  -5/ 10) 
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Pi^ure  19:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  the 
statistical  bias  variations. 

(wind  speed  =  0  kt;  cloudiness~5/lO) 
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Figure  20:   Total  temperature  variation  curves 
with  variation  in  starting  time  of  the  period. 
(wind  speed  =  0  kt;  cloudiness  =  C) 
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Figure  21:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  the 
statistical  bias  variations. 

(wind  speed  =  0  kt;  cloudiness  =  0) 
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Figure  22:   Total  temperature  variation  curves 

with  variation  in  starting  time  of  the  period. 

(2  kt^  wind  speed  =^15  kt;  cloudiness-^  5/10) 
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Figure  23:   Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  the 
statistical  bias  variations. 

(2  kt-wind  speed  ^.15  let;  cloudiness  «£  5/10) 
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Figure  24:   Total  temperature  variation  curves 
with  variation  in  starting  tirae  of  the  period. 
(wind  speed -A-.15  let;  cloudiness  *£10/l6) 
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Figure  25:  Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  statis- 
tical variations. 

(wind  speed .4-  15  let;  cloudiness  £-  10/10) 
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8.    48-hour  temperature  variation 

As  previously  mentioned 9   the  purpose  of  the  48-hour  tem- 
perature variation  analysis  is  to  examine  whether  the  statis- 
tical bias  occurs  for  periods  other  than  24  hours.   The  results 
of  this  analysis  are  presented  graphically  in  figures  26 
through  30. 

Figure  27  clearly  shows  an  effect  of  statistical  bias.  For 
this  figure  the  selection  criteria  were  windy  and/or  cloudy 
periods,  which,  according  to  Hisdals  should  have  their  highest 
temperatures  near  the  middle  of  the  periodo   This  is  what  ac- 
tually occurred,  but  only  five  48-hour  periods  met  the  speci- 
fied selection  criteria. 

In  figure  29 9  a  total  of  39  periods  met  the  criteria^  but 
the  actual  results  are  inconsistent  with  those  expected,  Accor= 
ding  to  the  discussion  of  section  1,  for  days  with  wind  speed 
less  than  or  equal  to  15  knots  and  cloudiness  less  than  or 
equal  to  five->tenths  the  temperature  trend  should  be  concave 
in  shape.   In  figure  29?  tn<§  statistical  bias  is  convex«   This 
departure  from  a  priori  expectation  may  be  due  to  the  random 
variations  included  in  the  total  temperature  eurve9  from  which 
the  bias  was  deduced*   This  idea  is  treated  more  thoroughly  in 
section  11 . 

The  last  figure  in  the  48-hour  set9  figure  30 B  illustrates 
the  mean  48-hour  temperature  variation  which  results  when  no 
particular  selection  criteria  are  specified*.   In  theory  these 
curves  should  have  negligible  statistical  biaSj,  considering  the 
discussion  in  section  1 a   This  theory  is  essentially  borne  out 
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in  figure  30,  where  the  mean  meteorological  variation,  obtained 
by  the  column  averaging  method,  is  very  nearly  the  same  as  the 
total  variation.   Thus,  upon  subtraction  of  the  meteorological 
from  the  total  variation,  the  resulting  statistical  bias  does 
not  appear  significantly  different  from  zero. 

It  should  be  noted  that  the  conclusion  of  a  zero  bias 
through  such  a  procedure  is  complicated  by  end  effects  due  to 
different  starting  times.   For  example,  the  first  two  three- 
hourly  temperatures  used  in  computing  the  0000  GOT  to  0000  GOT 
mean  curve  are  not  included  in  the  0600  GOT  to  0600  GOT  curve. 
Likewise,  the  last  two  temperatures  in  the  0600  GOT  to  0600  GOT 
curve  are  not  Included  in  the  0000  GOT  to  0000  GOT  curve. 
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Figure  26:  48-hour  temperature  variation 
curves  with  variation  in  starting  time 
of  the  period. 

(wind  speed ^15  kt;  cloudiness  ^5/10) 
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Figure  27:  Comparison  of  the  total  48-hour 
temperature  variation  with  the  meteorological 
and  statistical  bias  variations • 

(wind  speed  ="15  kt;  cloudiness^  5/10) 
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Figure  28:  48-hour  temperature  variation  curves 
with  variation  in  starting  time  of  the  period, 
(wind  speeds  15  let;  cloudiness^  5/10) 
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Figure  29:  Comparison  of  the  total  48-hour 
temperature  variation  with  the  meteorological 
and  statistical  bias  variation. 

(wind  speed*^15  let;  cloudiness  ^5/10) 
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Figure  30:  48-Hour  temperature  variation 
curves  with  variation  in  starting  time 
of  period(for  all  observations,  regard- 
less of  windspeed  and  cloud  amount). 
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9.    1 8-hour  and  72-hour  temperature  variations 

As  with  the  48-hour  temperature  variation,  the  purpose 
for  studying  the  18-hour  and  72-hour  variations  is  to  demon- 
strate that  a  statistical  bias  effect  is  not  limited  to  a 
24-hour  period. 

The  13-hour  temperature  variation  curves,  figures  31  and 
32,  clearly  depict  the  existence  of  a  statistical  bias.   In 
principle,  these  curves  should  be  concave,  since  the  criteria 
specified  are  wind  speed  less  than  or  equal  to  15  knots  and 
cloudiness  less  than  or  equal  to  five-tenths.   The  actual 
curves  are,  in  fact,  concave. 

In  figure  32  the  meteorological  temperature  variation  and 
the  statistical  bias  have  been  computed  by  means  of  the  two 
schemes  proposed  in  section  4;  that  is,  by  both  the  column- 
averaging  and  diagonal-averaging  processes.   The  resulting 
curves  are  quite  similar. 

Figures  33  and  34  depict  the  72-hour  temperature  varia- 
tion curves.   These  curves  are  designed  to  test  the  hypothesis 
that  the  statistical  bias  effect  exists  for  considerably 
longer  periods  than  24  hours.   However^,  they  do  not  support 
this  hypothesis.   In  fact9  there  is  evidence  of  an  upper  limit 
for  the  length  of  the  period  yielding  a  bias  effect.   A  possi- 
ble explanation  for  this  limitation  is  that  the  physical 
mechanism  causing  the  type  of  day  to  persist  for  a  long  period 
remains  long  enough  so  that  the  observations  on  either  side  of 
the  selected  period  also  meet  the  criteria.   If  this  were  the 
case,  no  statistical  bias  would  be  evidenced. 
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A  further  investigation  is  necessary  in  order  to  deter- 
mine which  synoptic  situations  at  Antarctica  are  associated 
with  long  sequences  of  similar  days. 

One  difficulty  which  arises  in  studying  a  period  as  long 
as  72  hours  is  that  the  number  of  periods  meeting  the  criteria 
are  relatively  few.   Thus,  the  statistical  significance  is 
greatly  reduced  compared  to  that  of  shorter  periods. 
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Figure  31:  1 8-hour  temperature  variation 
curves  with  variation  in  starting  time 
of  the  period. 

(wind  speeds! 5  let;  cloudiness^5/lO) 
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Figure  32:  Comparison  of  the  total  18-hour 
temperature  variation  with  the  meteorological 
and  statistical  bias  variations, 

(wind  speed -^15  ^t;  cloudiness  s£5/l0) 
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Figure  34:  Comparison  of  the  total  72-hour 
temperature  variation  with  the  meteorological 
and  statistical  bias  variations, 

(wind  speeds  15  kt;  cloudiness^  5/10) 
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10.   The  diurnal  temperature  variation,  based  upon 
all  observations  1956  -  1961 

Figure  35  is  a  chart  of  mean  diurnal  or  meteorological 
temperature  variation  curves,  based  on  all  three-hourly  sur- 
face observations  at  McMurdo  Sound,  during  the  dark  season, 
for  the  years  1956  through  1961.   The  first  six  of  these  curves 
are  annual  means,  based  on  observations  of  93  days.   The  last 
curve  is  the  six-year  mean  diurnal  variation,  based  upon  558 
days. 

Since  no  selection  criteria  were  specified  for  these  mean 
curves,  each  represents  a  diurnal  temperature  variation,  that 
is,  one  free  from  statistical  bias. 

There  are  two  pronounced  features  of  the  six-year  mean 
curve.   The  first  is  the  temperature  maximum  occurring  at  1500 
GOT.   This  maximum,  however,  may  have  no  particular  significance, 
since  it  is  caused,  in  part,  by  the  abnormally  high  1500  GOT 
mean  temperature  which  occurred  in  1956.   Moreover,  during  the 
years  1958  through  1961,  there  is  a  minimum  value  at  1500  GOT. 

The  second  feature  of  the  six-year  mean  diurnal  temperature 
variation  curve  is  the  maximum  occurring  at  0000  GOT.   This 
appears  to  be  a  regularly  occurring  phenomenon,  at  least  for 
the  years  included  in  this  study. 

In  order  to  discover  a  possible  explanation  for  the  0000 
GOT  maximum,  the  diurnal  variations  of  some  other  meteorolo- 
gical parameters  were  plotted  in  figures  36  through  42.   The 
first  of  these,  figure  36,  gives  an  estimate  of  the "mean  diur- 
nal variation  of  wind  direction."  In  order  to  simplify  the 
computation,  a  scalar-average  wind  direction  was  computed 
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instead  of  a  vector  average.   However ,  the  diurnal  course  of 
wind  direction  did  not  reveal  a  modal  value  at  0000  GOT.   On 
the  other  hand,  Hisdal  [1 ,  p. 172]  in  his  investigation  did, 
in  fact,  obtain  a  relationship  between  wind  direction  and  sur- 
face temperature,  using  a  histogram  of  wind  direction  in  con- 
Junction  with  the  diurnal  temperature  trend.   It  would  be 
interesting  to  consider  histograms  of  the  wind-direction  fluc- 
tuation associated  with  the  temperature  curves  of  this  study. 

The  next  parameter  examined  was  the  diurnal  variation  of 
cloudiness o   The  annual  means  and  six-yearly  mean  are  found 
in  figure  37.   This  set  of  curves  bears  a  very  close  resem- 
blance to  the  mean  curves  of  diurnal  temperature  variation, 
even  to  the  point  of  having  a  0000  GOT  maximum.   Once  again, 
the  conclusion  may  be  drawn  that  cloudiness  plays  a  very  im- 
portant role  in  the  mean  diurnal  temperature  fluctuation  during 
the  polar  night.  However,  it  cannot  be  concluded  that  the 
diurnal  fluctuation  of  cloudiness  is  the  cause  of  the  diurnal 
temperature  variation,  for  then  it  becomes  quite  difficult  to 
explain  why  there  should  be  a  systematic  periodicity  in  cloud- 
iness.  It  seems  more  reasonable  to  believe  that  the  two 
phenomena  are  due  to  some  common  cause 0 

A  relationship  between  cloudiness  and  temperature  would 
also  be  indicated,  if  a  comparison  were  made  between  figures  39 
and  19.   In  these  curves,  the  selection  criteria  were  calm  wind 
and  cloudiness  less  than  or  equal  to  five-tenths.   These  cri- 
teria were  chosen  in  order  to  minimize  temperature  changes 
due  to  turbulent  mixing  associated  with  higher  wind  speeds. 
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Comparing  the  mean  variations  in  figures  19  and  39 »  it  appears 
that  for  periods  with  calm  winds  a  maximum  or  minimum  of  tem- 
perature tends  to  lag  the  corresponding  cloudiness  by  three 
hours.   Table  7  summarizes  some  observations  of  this  lag  re- 
lationship. 

Table  7 

Time  (GOT)  of  occurrence  of  maxima/minima 
for  days  with  wind  speed  =  0  and  cloudiness 
less  than  or  equal  to  five-tenths 

Cloudiness  Temperature 

Max.     Min.  Max„      Min. 

2400  2100 

2100  1800 

1 800  1 500 

1500  1200 

1200  0900 

0900  0600 

0300  0300 

Since  neither  mean  wind  direction  nor  mean  cloudiness 
yielded  a  causative  explanation  of  the  diurnal  variation  of 
temperature,  another  parameter  was  investigated.   The  parameter 
selected  next  was  the  sea-level  pressure.   Its  annual  mean 
diurnal  variation  and  six-year  mean  are  found  in  figure  40. 
Also,  figure  42  shows  a  comparison  of  the  six-year  mean  diur- 
nal pressure  variation  with  that  of  the  corresponding  temper- 
ature variation,  as  well  as  the  variation  of  other  parameters. 

The  important  features  of  the  six-°year  mean  diurnal 
pressure  variation  is  that  it  has  a  maximum  value  at  1200  GOT 
and  a  secondary  maximum  at  0000  GOT,   This  result  agrees  with 
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the  results  of  a  recent  investigation  of  Antarctic  data  by- 
Carpenter  [10]  .   He  used  two  years  of  three-hourly  surface 
pressure  data  in  making  a  harmonic  analysis  to  determine  the 
phase  and  amplitude  of  the  diurnal  and  semidiurnal  pressure 
waves.   Carpenter's  results  for  McMurdo  Sound  are  listed  in 
table  8  which  compares  favorably  with  the  pressure  variation 
in  figure  42. 

Table  8 

Mean  harmonic  components  of  the  diurnal 
pressure  wave  at  McMurdo  Sound,  after 
Carpenter  [10] 

Ampl.(mb)   Mean  time  of  max. 

Diurnal  Oscillation  .001         2118  GCT 

Semidiurnal  Oscillation        .076        1040  GCT 

In  view  of  the  fact  that  there  is  a  definite  diurnal  and 

semidiurnal  pressure  fluctuation  at  McMurdo  Sound,  the  next 

problem  is  to  determine  what  Influence  they  have,  if  any,  on 

the  diurnal  temperature  variation. 

Since  the  semidiurnal  variation  has  its  maximum  value 

near  1200  GCT,  it  is  readily  apparent  that  this  compressional 

pressure  wave  can  not  account  for  a  temperature  minimum  which 

occurs  at  the  same  time.   The  Poisson  equation  shows  that  for 

an  adiabatic  pressure  increase,  a  temperature  increase  will 

result.   Moreover,  the  diurnal  pressure  wave  (maximum  at  2118 

GCT)  has  an  amplitude  of  only  .001  mb.   This  pressure  increase 

would  result  in  a  temperature  increase  by  adiabatic  compression 

of  the  order  .0002P. 

Thus,  it  is  seen  that  neither  the  diurnal  nor  semidiurnal 
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Pigure35*  Mean  diurnal  temperature 

variation  during  the  polar  night 
for  the  years  1956  through  1961. 
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Figure36:  Mean  diurnal  variation  of  wind 
direction  during  the  polar  night  for  the 
years  1956  through  1961,  based  upon  all 
observations, . regardless  of  wind  speed 
and  cloud  amount. 
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Figure37:  Mean  diurnal  variation  of  cloudiness 
during  the  polar  night  for  the  years  1956 
through  1961,  "based  upon  all  observations, 
regardless  of  wind  speed  and  cloud  amount. 
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Figure  38:  Total  cloudiness  variation 
curves  with  variation  in  starting  time 
of  the  period. 

(wind  speed  =  0  kt;  cloudiness^  5/10) 
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Figure  39:  Comparison  of  the  total 
cloudiness  variation  with  the  meteor- 
ological and  statistical  bias  variations. 

(wind  speed  =  0  kt;  cloudiness^-  5/1 0) 
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Pigure40:  Mean  diurnal  variation  of  atmospheric 
pressure  during  the  isolar  night  for  the  years 
1956  through  1961,  based  upon  all  observations* 
regardless  of  wind  speed  and  cloud  amount. 
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Figure  41:  Comparison  of  the  total 
atmospheric  pressure  variation  with 
the  meteorological  and  statistical 
bias  variations. 


(wind  speed  =  0  kt;  cloudiness  =  0) 
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Figure  42:  Diurnal  variation  of  temperature, 
atmospheric  pressure,  wind  direction,  and 
cloudiness  during  the  polar  night.  These 
curves  are  based  upon  six  winters  of 
observations,  1956  through  1961 • 
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pressure  variations  can  account  for  the  0000  GOT  maximum  in 
the  diurnal  temperature  variation.   The  diurnal  pressure  wave 
is  not  of  sufficient  magnitude  to  affect  the  diurnal  tempera- 
ture variation,  and  the  semidiurnal  pressure  wave  is  12  hours 
out  of  phase  with  the  temperature. 

With  the  foregoing  physical  factors  eliminated  as  rea- 
sonable explanations  for  the  0000  GOT  maximum  in  the  daily 
temperature  curves  considered  previously,  it  is  necessary  to 
again  consider  a  statistical  explanation. 

A  climatological  statistical  effect  has  been  proposed  as 

an  explanation  by  van  der  Bijl^who  suggested  that  the  seasonal 

temperature  trend  could  produce  such  a  maximum.  A  schematic 

version  of  the  seasonal  effect  is  shown  in  figure  43. 

Figure  43 

Schematic  illustration  of 
the  seasonal  trend  effect 
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A  24-hour  segment  of  the  seasonal  temperature  trend, 
denoted  AB,  is  enlarged  and  shown  as  a  dashed  line  in  the 
upper  right-hand  corner  of  figure  43»   The  solid  lines  in  the 
enlarged  portion  show  the  segment  as  corrected  by  the  Lamont 
procedure. 

If  the  entire  temperature  trend  for  the  cooling  season 
were  divided  into  24-hour  segments,  each  would  have  a  concave 
curvature,  because  the  over-all  trend  is  concave .  Now,  if 
these  24-hour  segments  were  averaged  for  the  entire  cooling 
season,  the  mean  would  also  be  concave.   This  concavity  would 
give  the  appearance  of  a  diurnal  temperature  fluctuation,  such 
as  that  observed  in  the  present  investigation,,         * 

Figure  35  shows  that  the  concavity  of  the  mean  diurnal  or 
meteorological  temperature  variation  produces  a  temperature 
range  of  0.5F.   It  is  of  interest  then  to  determine  whether 
the  seasonal  temperature  trend  actually  has  enough  curvature 
to  produce  0.5F  in  a  24-hour  period.   In  pursuing  this  inves- 
tigation, the  seasonal  temperature  trends  were  plotted  for 
each  of  the  years  1956  through  1961.  Also9  the  six-year  mean 
trend  was  computed.   These  curves  are  shown  as  figures  44 
through  50. 

Examination  of  the  seasonal  temperature  trends  reveals 
that  they  are  not  as  smooth  as  was  expected.   Alsos  the  cur- 
vature from  the  beginning  to  the  end  of  the  dark  season  is 
not  of  sufficient  magnitude  to  cause  a  five-tenths  degree 
temperature  range  in  a  24-hour  segment.   Therefore ,  a  temper- 
ature-variation model  which  would  increase  the  curvature  of 
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the  average  24-hour  segments  would  substantiate  van  der  Bljl's 
proposal  and  yet  produce,  upon  averaging,  the  long-term  mean 
seasonal  temperature  curve. 

This  curvature  increase  may  be  accomplished  by  considering 
that  actual  observations  of  the  seasonal  temperature  trend  in  a 
given  year  are  as  shown  in  figure  51 .   In  this  figure  the 
dashed  line  represents  the  long-term  seasonal  trend,  which  is 
the  same  as  in  figure  43.   The  solid  line  represents  the 
unsmoothed  temperature  trend,  which  shows  that  relatively  long 
cooling  periods  during  the  dark  season  are  interrupted  by  short 
periods  of  above  seasonal  temperature. 

The  effect  is  to  produce  a  series  of  "spiked"  temperature 
maxima  which  greatly  increase  the  mean  curvature  of  the 
average  24-hour  segments  during  the  prolonged  cooling  segments. 
With  this  refinement  the  mean  diurnal  temperature  variation 
may  have  enough  curvature  to  produce  a  range  of  0.5F. 

A  more  complete  study  of  the  seasonal  temperature  trend 
is  necessary,  however,  before  its  effect  on  the  curvature  can 
be  validated.   In  particular,  if  the  shape  of  the  summer  tem- 
perature trend  is  as  shown  in  figure  51 9   and  the  total  mean 
diurnal  temperature  variation  for  the  mid-summer  months  has  a 
minimum  at  0000  GOT,  then  the  seasonal  trend  effect  would  have 
some  merit  for  the  summer  months  also. 
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Figure  51 

Schematic  model  of  a  proposed  seasonal 
temperature  trend  for  the  polar  regions 
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11 .   Discussion  of  random  errors 

Throughout  this  study,  random  errors  have  been  mentioned 
as  an  Important  factor  in  the  erratic  behavior  of  various 
curves.   For  instance,  the  2100  GOT  maximum  in  the  statistical 
bias  curve  of  figure  7  was  considered  to  be  caused  by  the 
randomness  of  some  of  the  data  used  in  computing  the  mean 
values.   The  following  discussion  will  illustrate  that  the 
2100  GOT  was,  in  fact,  caused  by  random  errors.  Also,  a  com- 
putational scheme  will  be  suggested,  which  reduces  the  effect 
of  these  errors. 

The  notation  used  is  as  follows: 

M]  =  temperature  due  to  the  meteorological 

influence  at  observation  time  i  =  0,..,  n,..,4 
Sj  =  temperature  due  to  the  statistical  bias 
t^  =  total  temperature  =  Mn  +  St  +  e,^ 
m  =  0, ....,4 
n  =  0, ... ,3 
emh=  temperature  due  to  the  random  error 
C2  =  variance  of  the  single- observation  random 
error  =  emh  . 
Subscripts: 

m  =  row  in  matrix  ; 
n  =  column  in  matrix, 

1,3=  chronological  order  of  the  temperature 
observation  . 
Consider  table  9  to  be  a  matrix  of  total  temperature 
similar  to  table  3.   The  meteorological  effect  is  constant 
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for  any  particular  column,  and  the  statistical  bias  is  con- 
stant for  any  particular  diagonal  row„   The  errors  are  assumed 
to  be  random  and  uncorrelated  in  both  the  vertical  and  the 
diagonal  directions. 

Table  9 

Generalized  Theoretical  Table 
of  Mean  Total  Temperatures 
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Using  the  column-averaging  method,  as  explained  in  sec- 
tion 4,  the  following  totals  result  from  using  the  temperatures 
in  table  9: 

Column  1:   4Mq  +  1-  Sj  +  £ewo 


)-.0 

3 


m-o 

a. 


Column  2:     4  M,  +  ^>_  Sj    +  ' 

J=0  m^c 


mi 


m-a 


Column  3:    4M*  +  ^  5,   +     ]>JW 
4:    4M3+   £s;    ♦      £c„s. 


Column 

Next,  divide  by  four  to  obtain  the  column  mean  temper- 
atures: 

,3      -   3 
Column  1:   M0  t  4  X^j  4  *  le^ 


J-0 


m-o 


Column  2:       M,    +  i    zL5i   +    A     zLe 


»ti 
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Column  3:   M,  +  *  Z  Sj  ^k  Z^ 


4  Z—  ctvi  x 


Column  4:  Kl3  +  i  T5;  +  i-  ^f  e 


4  Z_   ^3 


Subtracting  the  column  averages  from  their  corresponding 
first-row  entry  yields  an  estimated  value  of  the  statistical 
bias  for  the  various  columns : 


Column  1:   (  M0  +  $0   +  eoo) 

j-0         mnzc, 


|5*  -  i(S,  +  Sl4S3)4  ^e„~4  (e,0  te^^ej 
Column  2:   (  M,  +  5,  +  e0,  ) 

_i il£ .-L~jz).=fl..  — / 

Column  3:   (  Mx  +  S^  +  e„.) 

f  ^  - +( So  -  S,  +  53)  4  i  e62  -  !fe,z  +  e,4 t32) 

Column  4:   (  M3   +  S3    +  e0?) 

-  (M34^:|:Sj  +1  f.em3) 

v     3    Tj-o ~   WnsO £___ 

I  S3  -  4-Cs«  +  s,  +sz)+  |e63  -4r(el3tcx34eJ3) 

The  next  step  will  be  to  determine  the  variance  of  the 
statistical  bias,  in  order  to  be  able  to  compare  it  with  the 
variance  of  the  statistical  bias  computed  by  di agonal - aver- 
aging. 

If  a" is  the  variance  of  a  single-observation  random  error, 
then  the  variance  of  the  first  column=average  equals 
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1       q-  Q  ^oo 
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The  variance  of  the  statistical  bias  is  then  computed  as 
shown  below.   Note  in  this  connection  that  Mq  and   ^  ;>  Sj 
are  considered  to  be  fixed  numbers  and  therefore  have  zero 
variance. 

For  column  1:  V«r  [(^CeJ'  4*(e»°  *'e"  f  e'°)j 

-  ft--* 

The  other  columns  treated  similarly  will  yield  the  same  re- 
sult, which  is  the  standard  deviation  of  the  S±a 

The  second  method  considered  in  section  4  is  that  of 
averaging  the  total  temperatures  diagonally,,  As  explained 
previously,  a  diagonal  average  will  produce  the  statistical 
bias,  without  having  to  subtract  the  average  from  the  total 
temperature.  As  will  be  proven  below^,  this  method  results 
in  a  more  accurate  estimate  of  the  statistical  bias  than  that 
arrived  at  through  a  column   average.   Using  general  notation 
we  have 


Diagonal   sum:         4  S;   +  Z  Mi    +    2_£;i 


_±  4 

Average:  S;   +    fi.M;  t±  7i.t\\ 

;--i  T     :^i 


r-l  '       ;--| 
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Variance   in  averages      rr 


(^V 


Note  that  the  diagonal  method  reduces  the  variance  in  the 
statistical  bias  by  a  factor  of  three ,  compared  with  the  var- 
iance resulting  from  the  column  method. 

An  important  point  that  must  be  emphasized  is  that  the 
reduction  factor  of  three  is  valid  only  when  the  temperature 
observations  are  six  hours  apart.   In  other  cases,  where  the 
observations  are  other  than  six  hours  apart ,  different  re- 
duction factors  are  valid, 

A  generalized  derivation,,  similar  to  the  one  for  the 

six-hour  spacing,  can  be  used  to  arrive  at  the  values  in 

table  10,  which  are  the  factors  by  which  the  variance  in  the 

statistical  bias  will  be  reduced  if  a  diagonal  average  is  used 

in  computing  the  bias,  rather  than  a  column  average » 

Table  10 

Reduction  factor ,  R9  for  the  error 
variance  in  the  statistical  bias 


At 

1 

OJ 

3 

4 

6 

12 

24/k 

R 

23 

11 

7 

5 

3 

1 

k-1 

At  =  time  between  observations 

-   (24/k)  hours 
R   =  Reduction  factor 

This  discussion  of  random  errors  has  shown  that  the  best 
method  for  obtaining  the  statistical  bias  is  by  averaging  the 
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temperatures  in  table  9  diagonally,,   A  similar  discussion 
would  also  prove  that  the  best  method  for  obtaining  the  meteor- 
ological influence  is  through  a  column  average. 

Figure  52  was  plotted  in  order  to  provide  a  graphical 
comparison  of  the  results  of  these  two  methods.   The  most 
striking  feature  of  the  figure  is  that  using  the  diagonal 
average  in  computing  the  statistical  bias  has  eliminated  the 
2100  GOT  maximum. 
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Figure  52:  Comparison  of  the  total  temperature 
variation  with  the  meteorological  and  statis- 
tical variations. 
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12.   Conclusions  and  recommendations 

The  purpose  of  this  investigation  was  to  determine  the 
nature  of  the  diurnal  temperature  variation  during  the  polar 
night  at  McMurdo  Sound,  Antarctica.   The  conclusion  based  upon 
data  for  the  years  1956  through  1 961  is  that  no  diurnal  tem- 
perature variation,  in  the  normal  sense,  exists  during  the 
polar  night.   The  maxima  that  appear  in  the  daily  temperature 
trend  are  due  to  synoptic  fluctuations.   The  concave  curvature 
between  temperature  peaks  appears  to  reflect  that  of  the  sea- 
sonal cooling  curve. 

It  is  recommended  that  the  column- averaging  procedure 
should  be  used  only  for  determining  the  meteorological  compon- 
ent of  the  diurnal  temperature  variation.   For  computing  the 
statistical  bias,  the  diagonal-averaging  method  should  be  used. 
With  this  scheme,  the  variance  in  the  statistical  bias  can  be 
minimized. 

For  further  study  it  is  suggested  that  investigations, 
similar  to  the  present  one,  be  made  for  the  following  situations 

(1)  summer  season,  when  the  sun  never  goes  below 
the  horizonj 

(2)  other  polar  stations,  both  in  the  winter  and 
summer  seasons. 
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